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Abstract: In this study we evaluate the influences of optical property 
assumptions on near-infrared diffuse correlation spectroscopy (DCS) flow 
index measurements. The optical properties, absorption coefficient (u a ) and 
reduced scattering coefficient (u s f ), are independently varied using liquid 
phantoms and measured concurrently with the flow index using a hybrid 
optical system combining a dual-wavelength DCS flow device with a 
commercial frequency-domain tissue-oximeter. DCS flow indices are 
calculated at two wavelengths (785 and 830 nm) using measured u a and u s ' 
or assumed constant u a and u s '. Inaccurate u s ' assumptions resulted in much 
greater flow index errors than inaccurate u a . Underestimated/overestimated 
u s ' from -35%/+175% lead to flow index errors of +110%/-80%, whereas 
underestimated/overestimated u a from -40%/+ 150% lead to -20%/+40%, 
regardless of the wavelengths used. Examination of a clinical study 
involving human head and neck tumors indicates up to +280% flow index 
errors resulted from inter-patient optical property variations. These findings 
suggest that studies involving significant u a and u s ' changes should 
concurrently measure flow index and optical properties for accurate 
extraction of blood flow information. 
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1. Introduction 

Near-infrared (NIR) light has been recently employed in the noninvasive acquisition of blood 
flow information from deep tissues (up to several centimeters), which is referred to as NIR 
Diffuse Correlation Spectroscopy (DCS) [1-4] or Diffuse Wave Spectroscopy (DWS) [5,6]. 
DCS measures relative change of tissue blood flow (rBF) which has been extensively 
validated in various tissues through comparisons with laser Doppler flowmetry (LDF) [7], 
Doppler ultrasound [8,9], power Doppler ultrasound [10,11], Xenon-CT [12], fluorescent 
microsphere measurement [13], arterial spin labeled magnetic resonance imaging (ASL-MRI) 
[14,15], and to literatures [1,4,16-18]. DCS also provides a blood flow index for comparisons 
of longitudinal measurements and inter-subject variations [11,13,19,20]. The probing depth of 
NIR DCS (several centimeters) is significantly larger than those (several millimeters) of 
similar optical modalities such as LDF [21-23], Doppler optical coherence tomography 
(DOCT) [24], photoacoustic tomography (PAT) [25], and optical micro-angiography 
(OMAG) [26] . DCS is primarily sensitive to micro vasculature rather than large blood vessels 
(e.g., Doppler ultrasound measurement), and does not require radiation exposure (e.g., PET, 
Xenon-CT). Systems based on DCS provide portability, allowing for bedside monitoring 
utilizing short acquisition time (varying from 6.5 ms to several seconds) without expensive 
instrumentation [17,27-29]. Due to these features, usages of DCS expand continuously into 
new applications in various deep organs/tissues such as muscle [15,28,30-33], tumor 
[10,11,19,20,29,34-36] and brain [4,5,7-9,12-14,16,17,27,37-41]. 

The use of NIR light for deep tissue measurements stems from the exploitation of a 
spectral region (650-950 nm) wherein light absorption of the biological tissue is relatively 
low. When using NIR spectroscopy (NIRS) to detect optical properties of deep tissues, a pair 
of source and detector fibers is usually placed along the tissue surface with a distance of a few 
centimeters. NIR light generated by a laser emits into tissues through the source fiber and is 
detected by a photodetector through the detector fiber. Photon migration in tissue is now 
known to be a diffusive process [2,7]. During this migration, photons encounter absorption 
and, more commonly, scattering events. The probabilities of these events are described by an 
absorption coefficient, |i a , and a reduced scattering coefficient, |i s \ also referred to as the 
optical properties, intrinsic to the probed tissue volume. The penetration depth of NIR light in 
biological tissues is approximately half of the source-detector separation. NIR DCS flow 
measurements are accomplished by monitoring speckle fluctuations of photons emitted at the 
tissue surface. In non-muscular tissues moving red blood cells (RBC's) inside vessels are 
primarily responsible for these fluctuations [4,5,7-14,16,17,19,20,27,29,34^2], but 
complications such as tissue shearing and motion artifacts can arise for muscular tissues 
[28,32]. Blood flow indices and rBF can be calculated from the changes in the speckle 
patterns. Ensuing calculations of blood flow using DCS measurements include a dependence 
on the optical properties (|i a and |i s ') and are thus potentially influenced by variations thereof 
(see the details in Section 2). 
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DCS is not inherently capable of measuring absolute values of absorption and scattering 
coefficients. Solutions to this issue have typically been approached using two general 
methods: optical property assumptions or optical property measurements with separate 
instrumentation. Some studies have chosen to use the values of |i a and |i s ' from the literature 
[8,32], respective to the tissue type (e.g., brain or muscle), or assumed a constant |i s ' while 
examining changes in (i a [13,30,40]. These assumptions may be susceptible to deviations in 
optical properties that occur transiently, longitudinally, due to subject differences or from 
differences in literatures. A few of the recent studies have employed hybrid instrumentation 
allowing for measurement of both sets of information to extract accurate blood flow 
[9,19,29,43]. However, a generalization of potential flow index errors due to the inaccurate 
estimation of optical properties has not been investigated for the DCS flow measurements. In 
addition to optical properties, another potential influence on DCS flow indices is determined 
by selection of the laser wavelength. 

Our lab has recently built a hybrid instrument capable of the simultaneous measurement of 
absolute |i a , |i s ' and flow indices at multiple wavelengths, through combining a commercial 
frequency-domain NIR tissue-oximeter, the Imagent (ISS, Inc., IL, USA) [44,45], and a 
custom-made NIR DCS flow-oximeter [31,33,41]. This newly developed hybrid instrument 
allows us to quantify the influences of optical properties on DCS flow indices measured at 
different wavelengths. In this study, homogeneous liquid phantoms with controlled variations 
of optical properties were created, attempting to isolate the influence of each optical property 
parameter (i.e., |i a or lx s ') on DCS flow indices. The usage of tissue-like phantoms for 
instrument calibration and experimental validation of NIRS and DCS techniques is common 
[2,4,44-51]. In DCS measurements, the dynamic scatterer motions (typically 
micro vasculature RBC's) are best modeled by Brownian diffusion as opposed to random 
ballistic flow, which has been determined empirically, but for reasons currently unknown 
[4,7,8,10,11,20,27,40]. An effective Brownian diffusion coefficient is calculated as the blood 
flow index when measuring in biological tissues and is usually distinct from the conventional 
Brownian diffusion coefficient predicted by Einstein [52]. However, when utilizing liquid 
phantoms with Intralipid particles to provide Brownian motion, the two diffusion coefficients 
are expected to be equivalent. Through this special case using liquid phantoms, DCS flow 
indices calculated using measured or assumed optical properties can be compared to the 
Einstein prediction (as a true flow index) for Brownian particles suspended in liquid. 
Measurement errors are then determined through these comparisons for DCS flow indices at 
different wavelengths. 

Simultaneous measurements of optical properties and blood flow indices are essential for 
extracting accurate hemodynamic information in tissues with transient, longitudinal and inter- 
subject differences in optical properties. To this end, we show a clinical study using the 
hybrid instrument to accurately quantify tissue optical properties and blood flow indices in 
head and neck tumors. The measurement errors in tumor blood flow indices induced by 
potential inaccurate estimations of tissue optical properties are ultimately discussed and 
compared to the phantom study results to determine the in vivo applicability thereof. 

2. Methods and materials 

2.7. Diffuse correlation spectroscopy (DCS) for blood flow measurement 

DCS flow indices are quantified by a dual -wavelength DCS system [31] with two long 
coherence length continuous-wave (CW) NIR laser sources at 785 and 830 nm (100 mW, 
Crystalaser, Inc., NV, USA). The DCS sources emit light alternately into the tissue via two 
multi-mode optical fibers bundled at the same location on the tissue surface (see Fig. la). 
Four detector fibers are tightly bundled and placed on the tissue surface at a distance of 1.5 
cm away from the source fibers, and each is connected to a single photon -counting avalanche 
photodiode (APD) (PerkinElmer Inc., Canada). The outputs of 4 APDs are sent to a 4 -channel 
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autocorrelator board (Correlator.com, NJ, USA) producing normalized light intensity 
temporal autocorrelation functions (g 2 ) which are averaged to improve the signal-noise-ratio. 
The averaged g 2 from DCS is related to the normalized electric field temporal autocorrelation 
function (g x ) through the following Siegert relation [53]: 

g 2 (r,T) = \ + f3\ gl {r,T)\ (1) 

where x is the delay time, f is the position vector, and P depends on laser stability and 
coherence length and the number of speckles detected. 




Fig. 1. (a) Hybrid fiber-optic probe comprised of: two DCS source fibers (785 and 830 nm) and 
four bundled DCS detector fibers separated by a distance of 1.5 cm; eight ISS source fibers 
(780 and 830 nm) and ISS detector fiber separated by distances of 2.0, 2.5, 3.0 and 3.5 cm (2 
source fibers per separation distance). Note that two additional DCS detector fibers are shown, 
at 2.4 and 2.8 cm separations, but were not connected for this study, (b) Cartesian coordinates 
oriented for DCS source and detector (left) and liquid phantom setup (right) including: -9.5 L 
glass aquarium (30.5 cm x 21.0 cm x 15.0 cm), hybrid optical probe and holder, and lab stand, 
(c) Typical correlation function fitting from a phantom experiment (ji a (830 nm) = 0.05 cm" 1 , 
|V (830 nm) =10 cm" 1 ) with gi derived from g 2 measurements (gi m ) using Eq. (1) and gi 
calculated (g Jc ) using Eq. (3). 

Scatterer motion is directly associated with the unnormalized electric field temporal 
autocorrelation function (GO which obeys a correlation diffusion equation, derived rigorously 
elsewhere [2,46] and defined as follows for homogeneous media using a CW source (steady 
state): 

DV 2 -v Ma -^VjU s k 2 0 a(Ar 2 (r))^ (r,r) = -vS(r) (2) 

where D = v I (3jU s ) is the photon diffusion coefficient, v is the speed of light in the medium, 

k 2 is the wavenumber, S(r) is the source light distribution, and^Ar 2 (r)^ is the mean-square 

displacement of scatterers in time x. The position vector, r , denotes a general vector from a 
source to a point of detection. Note that Gi is the unnormalized version of gi, i.e., 
g 1 (f,r) = G x (r,r)/G 1 (r,0) . Scatterer movement for particles undergoing Brownian motion 

results in ^Ar 2 (r)^ = 6D B r , where D B is the effective Brownian diffusion coefficient. An a 

(0-1) term is added to account for not all scatterers being dynamic and is defined as the ratio 
of moving scatterers to total scatterers. The combined term, aD B , is referred to as the blood 
flow index in biological tissues and is commonly used to calculate the relative blood flow 
(compared to the baseline flow index before physiological changes). In contrast to tissue 
samples where scatterers may be static (e.g., organelle, mitochondria) or dynamic (moving 
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RBC), all scatterers in liquid phantom solutions (see Section 2.4) are considered dynamic with 
a »1 and the flow index is thus reported as simply D B . The homogeneous CW solution to Eq. 
(2) for semi-infinite geometry is 



where p is the source-detector separation, S 0 is source intensity, 



K 2 (r) = 3 MaMs + M Xa(Ar 2 (T)) , r x = [p 2 + {z-z 0 ) 2 ^ , r 2 = [p 2 + (z + z 0 + 2z b ) 2 ^ , 
z 0 = 1/ ju s , z b =2(l + R eff )l3ju s (l-R eff ) , R eff = -lA40n~ 2 + 0.710w _1 +0.668 + 0.0636^ and 



n « 1.33 (for tissues and phantoms) [7,46,54,55]. The R e ff term accounts for the mismatch 
between the medium and the air indices of refraction with n being the ratio between them. 

For semi-infmite geometry, the collimating laser source at (0, 0, 0) and detector at (p, 0, 0) 
are placed on the tissue surface with z = 0 (see Fig. lb). This solution (Eq. (3)) involves an 
extrapolated zero boundary condition including an isotropic source at z = z 0 and negative 
isotropic imaging source at z = -(z 0 + 2z b ). The position vector, f , from Eq. (2) considers the 
point source at (0, 0, z 0 ) and the negative imaging source at (0, 0, -(z 0 + 2z b )). The 
superposition of solutions to these two sources using infinite geometry provides the resulting 
Eq. (3) where now the semi-infinite boundary is modeled by the scalar parameter, p. Further 
details can be found elsewhere [46,55]. 

Flow index calculations begin with using Eq. (1) to first determine the p. Using the g 2 data 
at earliest x and letting gi »1 (i.e., g 1 (p,0) = G 1 (p,0)/G l (/?,0) = 1 ) lead to 
J3 = g 2 (p,r « 0)-l . Using DCS measured g 2 (p,r) , calculated p and Eq. (1), g x (p,r) is 
calculated for all x. Equation (3) is then used with the unknown parameter aD B (flow index) to 
fit the gi derived from DCS measurements (see Fig. lc). For a complete frame of DCS data 
acquisitions at two wavelengths, two flow indices are obtained sequentially. 

2.2. Frequency -domain spatially resolved NIRS for tissue optical property measurement 

Quantification of absolute (i a and |i s ' is performed by a frequency-domain multi-distance 
spatially resolved spectroscopy, i.e., the Imagent. Two wavelengths (780 and 830 nm) of a 
four-wavelength (690, 750, 780 and 830 nm) Imagent system are chosen to match the DCS 
lasers available (785 and 830 nm). The Imagent emits sinusoidally modulated light into tissue 
through 8 optical fibers (4 per wavelength) placed at four pre -determined distances (2.0, 2.5, 
3.0, and 3.5 cm) from a detector fiber bundle connected to a photomultiplier tube (PMT) (see 
Fig. la). Source light is modulated at 110 MHz resulting in detected AC, DC and Phase (cp) 
information from multiple distances. A simplified solution based on semi -infinite geometry 
for the photon diffusion equation exposes linear relationships between cp, logarithmic AC or 
logarithmic DC and spatial distances [44,45]. From the fitting slopes (S A c, S DC , S 9 ) of the 
linear relationships |i a and |i s ' can be extracted at each wavelength. Different source -detector 
separations generally provide measurements at different depths/regions based on diffusion 
theory [3]. However, the depth/regional difference has minimal effect on measurement of 
optical properties of homogeneous phantoms. 

2.3. Brownian motion of spherical particles in liquid phantoms 

As mentioned earlier, when using liquid phantoms with Intralipid particles to provide 
Brownian motion, the effective Brownian diffusion coefficient (flow index) measured by DCS 
is expected to be equivalent to the conventional Brownian diffusion coefficient predicted by 
Einstein [52]. In the present study, DCS flow indices are compared to Einstein diffusion 



V S 0 ( exp(-ir(r)r 1 ) exp(-K(r)r 2 ) 




(3) 
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coefficients for the estimation of measurement errors. The Brownian diffusion coefficient as 
defined by the Einstein-Stokes formula for spherical particles in liquid is 



where k B is the Boltzmann constant, T is the phantom temperature, R is the radius of the 
spherical particles and r| is the viscosity [52]. In order to obtain the temperature and viscosity 
parameters, a temperature sensor (Physitemp, NJ, USA) is attached near the optical probe and 
viscosity is measured using a viscometer (Brookfield, MA, USA). Viscosity is reported in 
units of cP (centipoise), where 1 cP = 1 mPa-s (millipascal- second) = 0.001 kg-irfV 1 [56]. 
The radius of Intralipid particles is estimated as 196 nm (see Section 2.5). 

2.4. Liquid phantoms with varied optical properties 

Liquid phantoms are comprised of distilled water, India ink (Black India 44201, Higgins, MA, 
USA) and Intralipid (30%, Fresenius Kabi, Uppsala, Sweden). India ink is used to manipulate 
the absorption coefficient of the phantom, ji a (k), where k is the laser source wavelength. India 
ink is first diluted to a 10% solution with distilled water. The 10% ink solution (instead of 
pure ink) is used to create phantoms. Intralipid provides particle Brownian motion and control 
of the reduced scattering coefficient of the phantom, |i s ' (k). Setup of the liquid phantom is 
shown in Fig. lb. A hybrid fiber-optic probe is placed on the surface of the liquid phantom 
solution contained inside a glass aquarium. A custom-made probe holder attached to a lab 
stand holds the probe at the center of the solution to simulate a semi-infinite geometry. 

For creating phantoms with |i a variation, a list of |i a (k) and a constant |i s ' (k) are chosen 
(see details in Section 2.5). The |i a (k) and |i s ' (k) of distilled water, 10% ink solution, and 
30% Intralipid are first determined. These values in combination with titration equations 
provide the necessary volumes of water, ink and Intralipid to achieve desired phantom optical 
properties [57]. The subscripts "ink", "Intralipid", and "water" are used in this paper to 
denote 10% ink solution, 30% Intralipid, and distilled water, respectively. The |i a ink (k) is 
derived from spectrometer (Beckman Coulter, CA, USA) measurements. Since the 10% ink 
solution is out of the measurable range of the spectrometer, further dilution is performed to 
get a 0.025% ink solution for spectrometer absorbance measurements. Absorbance measured 
from multiple 0.025% ink solution samples are averaged and converted to the u. aink (k) [57]. 
The |i a wa ter (^) is taken from the literature [58]. The |i a intralipid (fy is assumed to be equivalent to 
that of water. Distilled water and 10% ink solution are both assumed to have no contributions 
to the phantom ji s ' (k), i.e., u- s ' wa ter (^) = m'ink (A) = 0 cm" 1 . The |i s ' (k) of 10% Intralipid is 
calculated using a Mie theory approximation [59]. The theory and details including the 
Intralipid particle radius and refractive index were described in the original derivation [59] 
which has been extensively used for quantification of Intralipid-based liquid phantoms 
[2,7,57,60,61]. For 30% Intralipid used in this study, ji s 'intraiipid (k) is the |i s ' (k) of 10% 
Intralipid multiplied by a factor of three [57]. 

2.5. Phantom experimental protocols and data analysis 

Ha variation. In this protocol, variation of ji a was induced while maintaining a constant |i s '. 
Thirteen steps were performed to cover ji a (830 nm) from 0.05 to 0.20 cm" 1 with a step size of 
0.0125 cm" 1 (i.e., |i a (830 nm) = 0.05, 0.0625, 0.075, 0.20 cm" 1 ) and \i s ' (830 nm) = 10 
cm" 1 . Prior to beginning the Imagent requires calibration to a phantom of known optical 
properties. During this process, corrections are made to account for the efficiency of optical 
coupling among the lasers/detector, optical fibers, and phantom [45,62]. The Imagent was 
calibrated to a liquid phantom of equivalent composition and optical properties at the 
midpoint (step 7) of the experimental range, i.e., ji a (830 nm) = 0.125 cm" 1 and |i s ' (830 nm) = 
10 cm" 1 . The combined probe was then placed upon a liquid phantom at the lowest optical 
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property step, i.e., |i a (830 nm) = 0.05 cm -1 . For each of 13 steps the following actions were 
taken. Desired volume of ink solution was added to the liquid phantom, mixed, and allowed to 
stabilize for 10 minutes. Viscosity data was attained during this time by extracting three 500 
ul samples. The three samples were carefully taken from the left, right and middle of the 
solution (at the surface) to minimize spatial variations without submerging the pipette. Room 
light was then turned off and the experimental setup was covered with black plastic to reduce 
ambient light. Measurements were taken by the hybrid optical instrument and temperature 
sensor for a 5 minute interval. 

H k variation. Variation of |i s ' immediately followed the performing of ji a variation 
experiments. Between experiments the phantom from ji a variation was disposed of and 
replaced with a new initial phantom for |i s ' variation. The hybrid probe was cleaned with 
alcohol pads and repositioned on the surface of the second phantom for ji s ' variation study. 
Neither Imagent nor DCS were shut down between protocols. Variation of |i s ' was induced 
while maintaining a constant ji a . A scattering range of |i s ' (830 nm) from 4 to 16 cm" 1 with a 
step size of 1 cm" 1 (i.e., u s ' (830 nm) = 4, 5, 6, . . ., 16 cm" 1 ) and ji a (830 nm) = 0.125 cm" 1 was 
performed over thirteen steps. The volume of Intralipid required to increase |i s ' as desired 
could potentially reduce |i a as well as influence the level of probe submersion. To remedy the 
first difficulty, ink was added with the Intralipid to maintain the ji a of the phantom. For the 
second, an equivalent amount of phantom solution was removed after being mixed at each 
step. Viscosity, temperature and hybrid optical measurements were taken similarly to those 
during ji a variation. 

Data analysis and presentation. Each 5 minute interval measurement involves post 
calculations of interval averages of |i a and |i s ' at each wavelength along with the temperature, 
sample average of viscosity, and three diffusion coefficients (D B 's). Data between intervals 
(i.e., adding/taking solution, stirring) are excluded from data analysis. The ji a and |i s ' are 
measured by Imagent using the spatially resolved (slope) method (see Section 2.2) and 
averaged over the 5 minute interval. The averaged ji a and u s ' are then used as known 
parameters in fitting DCS measured g 2 's using Eq. (1) and (3) to produce two distinct DCS 
D B 's, which are distinguished with subscripts describing the optical property values used for 
calculation of D B . The first DCS D B (D B -mid) * s calculated with the averaged u a and u s ' 
measured from the middle interval, i.e., ji a (830 nm) = 0.125 cm" 1 and |i s ' (830 nm) = 10 cm" 1 , 
serving as the diffusion coefficient ignorant of optical property variation. Using the constant 
ji a and u s ' from the middle interval results in overestimations of ji a and u s ' during early 
intervals and underestimations at later intervals, thus causing errors in calculation of flow 
indices. The second DCS D B (D B . dynamic ) is calculated using the averaged ji a and |i s ' measured 
from the corresponding interval, providing the best evaluations of DCS flow indices. These 
DCS D B calculations are repeated for both sets of wavelengths. The optical properties from 
the Imagent measurements at 830 and 780 nm are used in calculations of DCS D B s at 830 and 
785 nm, respectively. The influence of wavelength mismatch (5 nm) between 780 and 785 nm 
is considered to be minor. The third D B (D B . Einstein ) is calculated using Eq. (4) with the 
estimated particle radius, interval averaged temperature and three sample averaged viscosity. 
The estimated particle radius of 196 nm is determined to exhibit the least errors between the 
measured DCS flow indices (D B . dynamic ) and calculated D B . Einstein at the calibration point (step 
7). This estimation falls within the range of Intralipid particle size reported in the literatures 
[2,59]. 

Measurement errors are characterized by calculating percentage errors between the 
measured values and predictions. P- values from Student t-tests are included for comparisons 
of measurement errors and the criterion for significance is p < 0.05. Results are presented in 
figures and tables to visualize measurement variations, differences between expected and 
obtained values, and the optical property influences on DCS flow indices. Data are depicted 
as mean + standard deviation (SD) in tables and SD is illustrated by error bars in figures. 
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2.6. In vivo quantification of head and neck tumor hemodynamics 

The possible errors of assuming constant optical properties in calculation of flow indices may 
be more readily visualized through analysis of real tissue measurements. In order to evaluate 
such influences in in vivo measurements, tissue hemodynamic properties of head and neck 
tumors in 10 patients were measured using a hybrid optical instrument similar to that in the 
phantom study. Institutional review board (IRB) approval was given by the University of 
Kentucky and consent forms were obtained before subject participation. Only subjects with 
Stage III-IVb Squamous Cell Carcinoma of the Head and Neck (SCCHN) were included. 
Neck lymph nodes measuring more than 1 cm and clinically thought to be involved by tumor 
were selected to study. 

The Imagent used 690 and 830 nm source wavelengths whereas DCS used 785 and 854 
nm lasers. Other wavelengths used for phantom studies were not available for the tumor 
study. Thus, data from 830 nm for Imagent and 854 nm for DCS are analyzed for this tumor 
study as these wavelengths provide the best match. By contrast, the second pair of 
wavelengths (785 nm versus 690 nm) is excluded for data analysis due to the significance of 
wavelength mismatch (95 nm). The Imagent source -detector separations used (2.0, 2.5, 3.0, 
and 3.5 cm) are identical to the phantom study. DCS utilized 3 detector fibers at 1.5, 2.4, and 
2.8 cm separations from the two source fibers. The probe was held by hand and secured on the 
subjects in the center of the area identified as tumor node while tumor optical properties and 
DCS flow data were obtained for -2 minutes. DCS data from the 2.8 cm separation are 
examined, comparable to the tissue region/depth probed by the Imagent. Using different sets 
of optical properties measured by the Imagent, four DCS flow indices for each of 10 subjects 
are calculated and then averaged over the 2-minute measurement interval. The |i a and |i s ' are 
averaged over the measurement duration (2 minutes) for each subject and used in calculating 
the first DCS blood flow index (aD B -d y namic)> which is considered as a true flow index. The 
minimum, mean and maximum ji a and |i s ' over 10 subjects are determined and used for 
calculating the respective remaining three DCS blood flow indices (aD B -mi n , aD B - me an, anc * 
aD B -max respectively) for comparisons with the true flow index (aD B . dynamic ). Data in figures 
are presented by interval mean + SD, where SD is depicted by error bars. 

3. Results 

3.1. ju a variation 

In order to evaluate the influence of u a variation on flow indices, thirteen steps of liquid 
phantoms were performed to cover ji a (830 nm) from 0.05 to 0.20 cm" 1 with a step size of 
0.0125 cm" 1 while maintaining a constant |i s ' (830 nm) = 10 cm" 1 . For each step/interval of 
measurements over 5 minutes, the means and SDs of viscosity (from three samples), 
temperature and calculated D B . Einstein are displayed as data sets (means) and error bars (SDs) in 
Fig. 2. D B . Einstein (see Fig. 2c) is calculated using the measured temperature (see Fig. 2b) per 
interval along with the associated viscosity (see Fig. 2a) and estimated particle radius (196 
nm). 

The interval means and SDs (error bars) of u. a , |i s \ D B . Einstein , D B -mid, and D B . dynamic 
throughout the 13 steps of |i a variation are displayed in Fig. 3. The measured values of |i a (see 
Figs. 3a and 3d) and |i s ' (see Figs. 3b and 3e) at two wavelengths by the Imagent are 
compared to predictions calculated using spectrometer and Mie theory for the evaluation of 
measurement errors, respectively. Two DCS flow indices (D B . mid and D B . dynamic ) are compared 
to the D B . Einstein (as a true flow index) for both wavelengths (see Figs. 3c and 3f). The D B .mid or 
D B . dy namic at each wavelength is calculated using the DCS measurement with averaged ji a and 
|i s ' from the middle interval [ji a (830 nm) = 0.125 cm" 1 and |i s ' (830 nm) = 10 cm" 1 ] or from 
the corresponding interval. 
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Fig. 2. Stability of Brownian motion (D B -Einstein) during p, a variation, depicted as means + SDs 
(as error bars). Viscosity (a) is averaged from three samples corresponding to each step; 
temperature (b) and D B -Einstein (c) are averaged respectively over 5 -minute intervals. 
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Fig. 3. Measured optical properties with calculated D B -Einstein and DCS flow indices during [i a 
variation. Left panels represent measurements obtained by Imagent and DCS at 830 nm. Right 
panels show measurements from Imagent at 780 nm and DCS D B calculations corresponding to 
785 nm. All values shown, except predicted |i a (a and d) and |V (b and e), are interval means ± 
SDs. Left panels: (a) |i a (830 nm), (b) \i s ' (830 nm), and (c) D B . Einst ein (830 nm), D B . dynamic (830 
nm), D B . mid (830 nm). Right panels: (d) p, a (780 nm), (e) u/ (780 nm), and (f) D B . E i n stein (785 
nm), D B . dynamic (785 nm), D B . mid (785 nm). 

3.2. ju/ variation 

Similar to ja a variation, a scattering range of |i s ' (830 nm) from 4 to 16 cm" 1 with a step size of 
1 cm" 1 was performed over thirteen steps while maintaining a constant ji a (830 nm) = 0.125 
cm" 1 . Results for |i s ' variation are plotted in a similar fashion as |i a variation (see Section 3.1). 
The means and SDs (error bars) of viscosity (three samples), temperature and calculated DB- 
Einstein throughout |i s ' variation are shown in Fig. 4. The interval means and SDs of (i a , |i s ', D B . 
Einstein, D B -dynamic, and De-mid are displayed in Fig. 5. 
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Fig. 4. Stability of Brownian motion (D B -Einstein) during p, s ' variation, depicted as means ± SDs 
(as error bars). Viscosity (a) is averaged from three samples corresponding to each step; 
temperature (b) and D B -Einstein (c) are averaged respectively over 5 -minute intervals. 
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Fig. 5. Measured optical properties with calculated D B -Einstein and flow indices during p, s ' 
variation. Left panels represent measurements obtained by Imagent and DCS at 830 nm. Right 
panels show measurements from Imagent at 780 nm and DCS D B calculations corresponding to 
785 nm. All values shown, except predicted |i a (a and d) and |i s ' (b and e), are interval means ± 
SDs. Left panels: (a) |i a (830 nm), (b) \i s ' (830 nm), and (c) D B . Einst ein (830 nm), D B . dynamic (830 
nm), D B . mid (830 nm). Right panels: (d) p, a (780 nm), (e) u/ (780 nm), and (f) D B . E i n stein (785 
nm), D B . dynamic (785 nm), D B . mid (785 nm). 

3.3. Quantification of ju a and ju/ influences onflow indices 

Influence of and variations on D^-Einstein - Table 1 lists the means ± SDs and 
coefficients of variation (CVs) for viscosity, temperature, and D B -Einstein over the entire range 
of |i a and |i a variations, calculated based on the data shown in Figs. 2 and 4. The CV is 
defined as a percentage of SD/mean, indicating the variation of the mean values over steps. 
The CVs of temperature, viscosity, and D B -Einstein are less than 2.2%, indicating the minor 
influences of |i a and (i a variations on these variables. The D B -Einstein is thus used as a true flow 
index to evaluate the DCS flow measurement errors. 
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Table 1. Mean ± SD and CV of viscosity, temperature and D B .Einstein 



Variables 


|x a variation 




[i s ' variation 




Mean ± SD 


CV 


Mean ± SD 


CV 


Viscosity (cP) 


0.98 ± 0.02 


1.54% 


0.96 ± 0.02 


2.12% 


Temperature (°C) 


18.61+0.20 


1.05% 


18.80 + 0.12 


0.64% 


D B -Einstein (CITIES) 


1.05E-08 + 1.65E-10 


1.57% 


1.07E-08 + 2.25E-10 


2.11% 



Mean measurement errors in n a , and DCS flow indices. Table 2 lists the means ± 
SDs of absolute percentage errors in measurements of |i a , |i s \ and DCS flow indices (D B -d y namic 
and D B -mid) over the entire range of |i a and |i a variations, calculated based on the data shown in 
Figs. 3 and 5. Absolute percentage error is defined as (IEstimate-True|/True) X 100%. For u. a 
and |i s \ the Imagent measured values are considered estimates while the spectrometer and Mie 
theory, respectively, are used as true values. For flow indices, D B -d y namic and D B .mid are 
considered estimates and D B . Einstein as true. The measurement errors for ji a , ji s \ and D B . dy namic 
are small, averaging less than -7%, whereas those of D B .mid are found to be larger, averaging 
up to 12.89% and 49.63% for |i a and |i s ' variations, respectively. The influences of |i s ' 
variation can be seen to produce greater percentage errors on flow indices than those of |i a 
variation. In Table 2, mean measurement errors between wavelengths are also compared using 
2-sample unequal variance, two-tailed t-tests with significant differences considered for p- 
value < 0.05 and denoted with * prefix. Significant differences in mean measurement errors 
between wavelengths are found in |i a (p = 0.01) during (i a variation and in |i s ' (p = 0.04) during 
|i s ' variation. These differences between wavelengths are most likely associated with the 
intrinsic feature of the instrument (Imagent) in detection accuracy at different wavelengths. 
No significant differences in mean measurement errors between wavelengths are found for 
both D B . dyn amic and D B -mid. 



Table 2. Imagent/DCS measurement percentage errors at 780/785 nm (upper line) and 
830/830 nm (lower line) 







|i a variation (Absolute % Error) 


[is variation (Absolute 


% Error) 


Variables 


Mean ± SD 


p -value 


Mean + SD 


p -value 




780 nm 


3.39 + 3.07 


*0.01 


1.86 + 1.15 


0.23 




830 nm 


0.84 ± 0.95 


2.86 + 2.69 




780 nm 
830 nm 


1.93 + 1.23 
1.60 + 1.50 


0.54 


3.14 + 2.66 
5.29 + 2.27 


*0.04 


DB-dynamic 


785 nm 
830 nm 


5.52 + 3.69 
4.02 + 3.30 


0.28 


5.84+10.73 
6.58 + 6.16 


0.83 


De-mid 


785 nm 
830 nm 


12.89 + 12.00 
10.89 + 8.99 


0.64 


49.63 + 31.51 
46.76 + 27.44 


0.81 



*p-values < 0.05 

Table 3. The p-values for comparisons of the mean measurement errors between D B .d y namic 



and D B -mid 


DB-dynamic VS. De-mid 


|i a variation (p-value) 


[is variation (p-value) 


785 nm 


0.0525 


*0.0003 


830 nm 


*0.0205 


*0.0002 



*p-values < 0.05. 



Table 3 provides p-value results for comparisons of the mean measurement errors of D B _ 
dynamic and D B -mid at two wavelengths during |i a and |i s ' variations. It is apparent that there are 
significant (though it is borderline at 785 nm during ji a variation) differences between the 
measuring (D B .d yna mic) and assuming (D B -mid) variables; D B -d y namic is more accurate (with less 
measurement errors, see Table 2) than DB-mid- The much lower p-values for ji s ' variation as 
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compared with those for |i a variation again suggest |i s ' as a greater influence factor on DCS 
flow indices. 

Influence of n. ± and variations on DCS flow index. Visualization of the influence of 
optical property assumptions at both wavelengths is shown in Fig. 6, overlaying the results 
from both |i a and |i s ' variations over 13 steps (see Figs. 3 and 5). The percentage errors for |i a 
during ji a variation and for ji s ' during u. s ' variation are defined as [(u. a -mid - |i a -dynamic)/^ -dynamic] 
X 100% and [(|V-mid ~ M-s'-dynamicVM's' -dynamic] X 100%, respectively. The subscripts "mid" and 
"dynamic" correspond to assumed constant (middle -interval) and dynamic optical properties. 
For both variations, the percentage D B error between D B - m id and D B .d y namic for each interval is 
defined as [(D B -mid - D B . dynamic )/D B . dynamic ] X 100%. Larger estimation errors in optical 
properties (|i a and |i s ') generate larger percentage D B errors. Variations in |i s ' have a greater 
influence on percentage D B errors compared to variations in ji a . Trends in overestimation and 
underestimation of flow indices due to variations in |i a or |i s ' are different. Overestimating and 
underestimating |i a results in overestimating and underestimating flow indices, respectively, 
opposite of the trend for |i s '. Data for both wavelengths are in good agreement and show only 
minor differences. 
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Fig. 6. Inaccurate estimations (percentage errors) of |i a and |V result in corresponding 
percentage D B errors between D B -dynamic and D B -mid for both wavelengths. 

3.4. Influence of tissue optical properties on head and neck tumor blood flow index 

The means ± SDs of measured tumor optical properties (|i a and |i s ') and blood flow indices 
(aD B . dynamic , aDE.min, aD B _ mean , and aD B . max ) along with corresponding percentage errors for 10 
patients with head and neck tumors are shown in Fig. 7. Patients are shown in order of 
increasing aD B . dynamic (as true flow index), designated with a black line, for comparison of 
trend differences when using optical property assumptions (i.e., aDe.^, aD B . mean , and aD B . 
max). Note that the patient numbers represent indices to illustrate the trend rather than actual 
patient numbers corresponding to the measurement sequence. The mean optical properties 
over subjects are: |i a (830 nm) = 0.12 + 0.03 cm" 1 and (830 nm) = 7.80 + 2.64 cm" 1 . 
Maximum and minimum optical properties out of all subjects at 830 nm are indicated using 
the red and blue dots, respectively, in Figs. 7a and 7b. DCS blood flow indices calculated 
using DCS data at 854 nm with different optical properties are represented in Fig. 7c. Without 
considering the tissue optical property influence, the trends of flow indices (aDB.^, aD B _ mean , 
and aD B _ max ) are not the same as the true flow index (aD B . dynamic ). Percentage aD B errors are 
calculated between the aD B . dynamic (true) and the estimated aDB.min, aD B _ mean , and aD B _ max . 
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Large ranges of percentage errors are found for aD B estimates: aD B -min from -8.07 to 
278.15%, aD B . me an from -39.48 to 149.01%, and aD B . max from -70.26 to 22.59%. The 
tendency to overestimate or underestimate the blood flow indices follow the same trends as 
shown for |i s ' variation in Fig. 6, supporting that |i s ' has a greater influence on DCS flow 
indices than |i a . 
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Fig. 7. Tumor optical properties measured by Imagent at 830 nm and flow indices measured by 
DCS at 854 nm using 830 nm optical properties. Data were obtained from tumor region for 10 
subjects with head and neck tumors. All values, except percent errors (d), shown are means ± 
SDs. (a) \i a (830 nm), (b) \i s ' (830 nm), (c) aD B -dynamic (854 nm), aD B - min (854 nm), aD B - me an (854 
nm) and aD B - m ax (854 nm), and (d) % Error aD B . min (854 nm), % Error aD B . mea n (854 nm) and % 
Error aD B . max (854 nm). Patients are listed in order of increasing aD B . dyn amic (854 nm). 

4. Discussion and conclusions 

4.1. ju a and ju/ variation influences on D B . Einstein 

The Einstein- Stokes formula calculation, Eq. (4), provides the Einstein diffuse coefficient 
(DB-Einstein) for spherical particles moving in liquid phantoms. The D B -Einstein is determined by 
the temperature and viscosity of the liquid as well as the particle radius of Intralipid in the 
liquid phantoms. Only slight variations are exhibited in overall average temperature (CV < 
1.1%), viscosity (CV < 2.2%), and D B -Einstein (CV < 2.2%) during both |i a and |i s ' variations, as 
seen in Table 1 and Figs. 2 and 4. The particle radius of Intralipid should not change during 
both |i a and |i s ' variations. The changes in temperature are likely due to room temperature 
increase over the -4.5 hour experimental durations. These include contributions by heat from 
equipment in the confined room. The cause of small variations in viscosity is likely due to the 
measurement variations by the viscometer. With these small variations in temperature and 
viscosity, it is thus expected that the derived D B -Einstein from Eq. (4) is stable over the large 
variations of optical properties. 

On the other hand, increases of (i a (ink concentration) during |i a variation are expected to 
have no contribution to D B -Einstein, as ink provides no particle motion. Similarly, increases of 
|i s ' (Intralipid concentration) during |i s ' variation do not show significant influence on D B _ 
Einstein, which is expected as all scatterers (Intralipid particles) provide equivalent motion in 
liquid phantoms and the ratio of moving scatterers to all scatterers (a) remains unchanged (a = 
1). Due to the independence of optical properties and high stability throughout, D B . E i ns tein is 
considered reasonable as the true flow index for spherical particles moving in liquid 
phantoms. 
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4.2. Measurement errors of ju a , ju/, and D B . dynamic 

In agreement with expectations, ink contributes only to increasing the absorption of the 
phantom during ji a variation. Additions of equivalent amounts of ink per interval resulted in 
linear increases in ji a for both Imagent wavelengths (see Figs. 3a and 3d). Only minimal 
variations occurred in |i s ' at both wavelengths during |i a variation (see Figs. 3b and 3e). 
Calibration at the midpoint [|i a (830 nm) = 0.125 cm 4 , ji s ' (830 nm) = 10 cm" 1 ] may influence 
the variation patterns seen in Imagent measurements. Intralipid contributes linearly to |i s ' 
increases, as expected, during |i s ' variation for both Imagent wavelengths (see Figs. 5b and 
5e). The |i a stayed relatively constant with minimal variations during ji s ' variation (see Figs. 
5a and 5d). Measured |i a and |i s ' during both experiments are consistent with predictions from 
spectrometer measurements and Mie theory, respectively. More specifically, the |i s ' 
measurement errors were less than 6% (see Table 2) for both experiments and wavelengths, 
which are comparable to those obtained from the literature using the Mie theory estimation 
(see Section 2.4) [59]. Overall, average measurement errors of ji a and |i s ' were small during |i a 
(< 4%) and |i s ' (< 6%) variations (see Table 2), which are in agreement with those found in 
previous studies using frequency-domain spatially resolved NIRS [44,45]. For ji a variation, 
significant difference was found between measurement errors of |i a at 780 and 830 nm. For |i s ' 
variation, there was significant difference between measurement errors of |i s ' at 780 and 830 
nm. These differences are likely attributable to the intrinsic instrument (Imagent) feature in 
detection accuracy at separate wavelengths. 

Average measurement errors for D B -d yn amic (< 7%) compared to D B .Einstein at both 
wavelengths are similar to those obtained for optical properties (< 6%) during both 
experiments (see Table 2), suggesting the influence of optical properties on D B . 

4.3. Resulting D B errors from optical property assumptions 

When using assumed constant optical properties (i.e., middle -interval |i a and |i s ') to calculate 
DB-mid> mean D B .mid measurement errors during |i a (-13%) and |i s ' (-50%) variations (see 
Tables 2 and 3) were significantly higher than those of D B . dynamic (-7%). It is evident that the 
influence of |i s ' on the DCS flow index is much greater than that of |i a . This result is further 
supported by the great difference in p-values, where p-values during ji s ' variation are much 
less than p-values during ji a variation (see Table 3). Also, looking at Fig. 6, the range of D B 
percentage errors for inaccurate estimations of |i s ' is much wider than that for |i a . This result is 
expected due to DCS flow indices being derived from light speckle fluctuations, originated 
from photon phase shifts by dynamic scatterers. Upon examination of the K 2 definition (see 
Eq. (3)), |i s ' should have a more significant influence than |i a given the |i s ' 2 term along with the 
much larger scattering over absorption (i.e., |i s ' » |i a ) in biological tissues and the liquid 
phantoms. No significant difference was found between wavelengths in D B . dynamic and D B -mid 
measurement errors during both experiments (see Table 2). This indicates that wavelength 
may not be a critical factor in determining the importance of optical property influence on D B 
measurement, although further investigations using a large range of wavelengths are needed 
for making a solid conclusion. The trends of D B estimation errors when using D B -mid were 
found to be different between the |i a and |i s ' variations (see Figs. 3c, 3f, 5c, 5f and 6). For |i a 
variation, overestimated or underestimated |i a results in overestimated or underestimated D B . 
By contrast, for |i s ' variation overestimated or underestimated |i s ' results in underestimated or 
overestimated D B . Extreme examples of incorrect estimations of D B can be seen in Fig. 6. 
Overestimation errors of |i a up to -+150% during |i a variation resulted in percentage errors up 
to -+40% and underestimation errors up of — 40% resulted in percentage errors up to — 20%. 
When overestimation errors of |i s ' reach up to ~+175% during |i s ' variation, D B percentage 
errors were up to — 80%. For underestimation errors of |i s ' up to — 35%, D B percentage errors 
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reach up to — hi 10%. Note that these estimation errors in optical properties and resulting D B 
may be affected by the selection of phantom properties for calibration. 

4.4. In vivo tumor study data in comparison to phantom study results 

In the tumor study, measured |i a and |i s ' show large variations between subjects (see Figs. 7a 
and 7b). The range of variations, |i a (830 nm) from 0.07 to 0.16 cm" 1 and fi s ' (830 nm) from 
5.35 to 13.1 cm" 1 , is within the range studied using liquid phantoms. The influence of the |i s ' 
variations on flow indices was found to be greater than that of (i a , supporting the phantom 
study results. This is exemplified by the trends shown in Fig. 7c. The overestimation of 
optical properties (using maximum ji a and |i s ') leads to underestimation of DCS flow index 
(aD B - max ) and underestimation (using minimum ji a and |i s ') leads to overestimation of DCS 
flow index (aD B -min)- These are in agreement with the trends of D B estimation errors using 
inaccurate |i s ' in liquid phantoms (see Figs. 5c, 5f and 6). Percentage aD B errors range greatly, 
from — 70% up to ~ +280%, depending on optical properties assumed. Errors in flow indices 
(see Fig. 7d) produce an incorrect observation of trends in the aD B magnitudes among patients 
(see Fig. 7c). It is evident that lack of consideration for optical property influences can lead to 
invalid results in similar studies. 

4.5. Conclusions 

The advent of DCS technology as a safe and quick alternative for measurement of blood flow 
in deep tissues has brought the need to further investigate potential errors, notably by the 
assumption of constant optical properties, |i a and |i s '. The flow index produced by DCS 
measurement is based on a solution to the correlation diffusion equation which includes 
parameters of |i a and |i s '. Utilizing a novel hybrid optical equipment setup, capable of 
measuring all three parameters of interest (i.e., flow index, ji a , and |i s '), with liquid phantom 
experimental protocols has made it possible to perform this investigation. The present study 
evaluates the influences of tissue optical properties on DCS flow indices through isolated 
variations of |i a and |i s ' in liquid phantoms. It is found that the particle motions in liquid 
phantoms are not influenced by the variations in optical properties, and the usage of Einstein 
particle Brownian motion coefficient (D B . Einstein ) as true flow index is reasonable for 
comparison with DCS flow indices. During ji a and |i s ' variations, |i s ' has a much greater 
influence on DCS flow indices than |i a , regardless of the wavelengths used. Studies involving 
significant |i a and |i s ' changes should concurrently measure flow index and optical properties 
for accurate extraction of blood flow information in tissue. The flow index errors resulted 
from the optical property assumptions in the tumor study elicit such need for concurrent 
monitoring of optical properties. Incorporation of laser sources at wavelengths beyond those 
tested in this study may be the subject of future investigation. The range of optical properties 
tested in the phantoms may also be extended to encompass a wider variety of tissues. 
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